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Abstract 



The behaviour of the broken CP symmetry at finite temperature is ex- 

■ amined. This is achieved through the investigation of the neutron electric 
dipole moment dn induced by ^-term. By using thermal QCD sum rules, we 
find that below the critical temperature, the ratio | ^ | slightly decreases 

' but survives at temperature effects. This evolution implies that CP remains 

■ broken at finite temperature as required by Baryogenesis [1]. 

1 Introduction 

Recently the finite temperature behaviour of symmetries has gained considerable 
^ ' interest. The question of symmetry restoration is a non trivial phenomenon, since 
' it has been shown in [23, 24] that more heat does not necessarily imply more 
symmetry. Besides, the breaking of the symmetries has a profound implications in 
a particle physics and cosmology. The CP symmetry is certainly one of the most 
fundamental symmetries in nature. Besides its role in solving domain wall problem 
[25], it is a crucial ingredient to understand the matter-antimatter asymmetry [26]. 

In the three generation Standard Model, CP violation originates from the more 
obscure sector of the SM: the scalar part. It is parameterized, in the electroweak 
sector, by a single phase occuring in the Cabbibo-Kobayashi-Maskawa (CKM) 
quark mixing matrix [2] . CP violation could also originate from additional CP-odd 
four dimensional operator embedded in the following topological term "6'-term" in 
the QCD Lagrangian: 

Le = e^G,,G^\ (1) 
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which breaks P, T and CP. is the gluonic field strength, denotes its 
dual and as is the strong coupling constant. The GuyG^" quantity is a total deriva- 
tive which contributes to the physical observables only through non perturbative 
effects, induced by instantons. A non zero value of 6 may generate, in particular, 
a sizable neutron electric dipole moment (NEDM) which is related to the ^-angle 
by the following equation obtained within the framework of Chiral perturbation 
theory: 

, e ,m^ ^ 2.7 X IQ-^^^ [3] 
^ ^ U.2 X 10--^ [4] 

High precision Experiments have constrained the NEDM to (i„ < 1.1 x 10~^^ ecm [5], 
providing a stringent upper limit to 6' < 2 x 10"^° [6]. The difficulty to explain the 
smallness of 9 in the standard model is usually known as the "strong CP problem" . 
In this regard, several scenarios were suggested. The most elegant explanation is 
due to Peccei and Quinn [7], who identified 6 to the axion, a very light pseudo 
scalar boson arising from the spontaneous breaking of a global C/a(1) symmetry. 
This particle may well be important to explain dark matter puzzle providing a 
peace of information on the missing mass of the universe [8] . 

Our aim in this work is to investigate the behaviour of the CP symmetry break- 
ing at finite temperature and the thermal effects on the restoration of the strong CP 
problem. This is motivated by the possibility to restore some broken symmetries 
by increasing the temperature. 

In section 2, we perform the the calculations of the the 9 induced NEDM using 
thermal QCD sum rules. Section 3 is devoted to the discussion and qualitative 
analysis of the thermal effects on the CP symmetry restoration. 



2 NEDM from thermal QCD sum rules 

In order to derive the NEDM through the QCD sum rules techniques [9, 10, 11], 
we consider a Lagrangian containing the following P and CP violating operators: 

Lp,cp = -9,m^ E Qfi^mi + e^G.^G'^r (3) 

9q and 9 are respectively two angles coming from the chiral and the topological 
terms while is the quark reduced mass given by m»= ^"^^ . the physical phase 
is ^ = ^ + We usually start from the two points correlators in QCD background 
with a non-vanishing 9 in the presence of a constant external electomagnetic field 

n(g2) = I j clW^^ < 0|r{r/(a;)?7(0)}|0 >e,F ■ (4) 
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where rj{x) is the neutron interpolating current [12]: 

V = 2eabc{{(i^C^5Ub)dc + p{(i^Cub)-f5dc}, (5) 

and /? is a mixing parameter. To select the appropriate Lorentz structure, n(g^) is 
expanded in terms of the electromagnetic charge as: 

n(g2) = n(0)(g2) + en(i)(g2, F'^-') + 0{e^). (6) 

The first term Il^'^^q'^) is the nucleon propagator which includes only the CP-even 
parameters, while the second term 11^^^ (q^ , F^") is the polarization tensor which 
may be expanded through Wilson OPE as: J2Cn < 0\qTq\0 >6i,f, where P is an 
arbitrary Lorentz structure and C„ are the Wilson coefficient functions calculable 
in perturbation theory [13, 14]. Prom this expansion, we keep only the CP-odd 
contribution part. The electromagnetic dependence of these matrix elements is 
determined in terms of the magnetic susceptibilities /c, x and ^ [14], defined as: 

< Olqaf'^qlO >f = XegF''" < 0\qq\0 > (7) 
g < OlqC^qlO >f = KegF^"" < 0|gg|0 > (8) 
2g < 0|gG"^''g|0 >f = ^e^F'^'^ < 0\qq\0 > (9) 

Moreover, the 9 dependence of < 0|gPg|0 >$ matrix elements may be traced by 
considering the anomalous axial current [10] : 

THg < 0\qTq\Q >e= imj < 0|gPg|0 > +0(m^) (10) 

where the correction 0{mg) is negligible since >> m^r. 

Putting altogether the above ingredients and after a straightforward calculation 
[11], the following expression of U.^^\q'^ , F^'^) for the neutron is derived: 



n(-a = -^<0\qq\0>{Fa,q}[x{P+inAe,-e^)\n{—^) 
-4(/?-lfe,(l + i(2/. + 0)(ln(#)-l) ^ 



4' Vfij' -q^ 

-^((4/3^ - 4/3 + 2)e, + (3/3' + 2(3+ l)e„)^...], (11) 
z q 

with q = qn'-f^. 

The QCD expression (11) is confronted to the phenomenological parameterization 
YiPhen^_q2^ Written in terms of the Neutron hadronic properties. The latter is 
given by: 

nP.e„(_^2) ^ {Fa,q}{j4^^^ + ^^A^ + ...), (12) 

(r - m2)2 (g2 _ ^2 j 
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where m„ is the neutron mass, is the quark charge, the parameters A and 
A^, which originate from the phenomenological side of the sum rule, represent 
respectively a constant of dimension 2 and the neutron coupling constant to the 
interpolating current r){x). This coupling is defined via a spinor v as < 0\'r]{x)\n >— 

In the framework of QCD sum rules, the correlators at finite temperature are 

expressed in terms of the thermal Gibbs average of Wilson operator expansion 
[15, 18]. At relatively low temperature, where the system can be regarded as a non 
interacting gas of bosons, the thermal dependence of the vacuum condensates can 
be written as : 

< O' >T=< 0'>+l < 7r(p)|0>(p) > nnip (13) 

where e = ^p'^ + m^, ub = is the Bose-Einstein distribution and < > is the 
standard vacuum condensate (i.e. at T=0). In this approximation, we only kept 
the pion contributions, since in the low temperature region, the effects of heavier 
resonances (F = rj, ..etc) are dumped by their distibution functions ~ [17]. 
To compute the pion matrix elements, we apply the soft pion theorem given by: 

1 777^ 

< n{p)\Oy{p) >= -- < 0\[QI [Qt,OW > +0(-rf ), (14) 
where A is a hadron scale and is the isovector axial charge defined by: 

Qt^ J d?xq{x)-iQ-i5—q{x). (15) 

Direct application of the above formula to the quark and gluon condensates shows 
the following features [16, 17]: 

(i) Only < qq> IS sensitive to temperature. Its behaviour at finite T is given by: 

<gg>rc^(l-^)<gg>, (16) 



where v?(T) = j^B{^), B{z) = ^ dy ^^y_^ - and is the pion decay constant 
(/tt ~ 93Mey). The variation with temperature of the quark condensate < qq >t 
results in two different asymptotic evolutions, namely: 

< qq >tC^ (1 - — ) <qq> 

for < 1, 



nnipi rp2 

< qq >T~ (1 - ^ e t )<qq> 
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for > 1. 

(ii) The gluon condensate is nearly constant at low temperature and a T dependence 
occurs only at order T^. 

The determination of the ratio ~ sum rules at non zero temperature is now 
easily performed by applying Borel operator to both parameterizations of the Neu- 
tron correlation function shown in (11) and (12). Then finite temperature effects 
are introduced via the procedure discussed above. Finally, by invoking the quark- 
hadron duality, we deduce the final sum rules of the 6 induced NEDM at finite 
temperature: 

(17) 

where M represents the Borel parameter. 

In order to get rid of the infrared divergence, the value of (3 has been set to 1 in 
(17). The Thermal evolution of the coupling constant and the mass of the neutron 
were determined from the thermal nucleon sum rules [17]. 

Within the dilute pion gas approximation, Eletsky has shown that the contri- 
bution induced the pion- nucleon scattering has to be considered [19]. It enters the 
nucleon sum rules through the coupling constant q-jtmn, whose values lie within the 
range 13.5-14.3 [20]. 

Numerical analysis is performed with the following input parameters: the Borel 
mass has been chosen within the values = 0.55 — O.TGeV"^ which correspond to 
the optimal range (Borel window) in the ~ sum rule at T = [11]. For the x 
^ susceptibilities we take x = -5-7 ± 0.6GeV-^ [22] and { = -0.74 ± 0.2 [?]. As 
to the vacuum quark condensate appearing in (17), we use its standard values [9]. 

3 Analysis and Conclusion 

We have established the relation between the NEDM and 9 angle at non zero 
temperature from QCD sum rules. We find that the behaviour of the ratio^ is 
connected to the thermal evolution of the pion parameters f^, and of Qt^nn- 
By analyzing the ratio as a function of T in the region of validity of thermal sum- 
rules [0,Tc], we learn that | ~- \ decreases smoothly with T (about 16% variation 
for temperature values up to 200 MeV) but survives at finite temperature. This 
means that either the NEDM value decreases or 9 increases. Consequently, for a 
fixed value of 9 the NEDM decreases but it does not exhibit any critical behaviour. 
Furthermore, if we start from a non vanishing 9 value at T = 0, it is not possible 
to remove it at finite temperature. We also note that | ~ \ grows as or x 
susceptibility increases. It also grows with quark condensate rising. However this 
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ratio is insensitive to both the ^ susceptibihty and the couphng constant g-nNN- 
We notice that for high temperatures, the analysis of | ^ |= /(^) exhibits a 
brutal increase justified by the fact that for T beyond the critical value Tc, at 
which the chiral symmetry is restored, the constants /„■ and q-^nn become zero 
and consequently the ratio ^ behaves as a non vanishing constant. The large 
discrepancy between the values of the ratio for T < T^. and T > T^. may originate 
from the other contributions to the the spectral function which have been neglected, 
such as the scattering process A'" + tt — > A. These contributions, which are of the 
order T^, are negligible in the low temperature region but become substantial 
for T > Tc- Moreover, this difference may also be due to the use of soft pion 
approximation which is valid essentially for low T (T < T^). Therefore it is clear 
from this qualitative analysis, which is based on the soft pion approximation, that 
temperature does not play a fundamental role in the suppression of the undesired 
^-term and hence the broken CP symmetry is not restored [1] . Indeed, some exact 
symmetries can be broken by increasing temperature [23, 24]. The symmetry non 
restoration phenomenon, which means that a broken symmetry at T=0 remains 
broken even at high temperature, is essential for discrete symmetries, CP symmetry 
in particular. Indeed, the symmetry non restoration is a crucial ingredient in 
solving the domain wall problem [25] and to create the baryon asymmetry in the 
early universe (BAU) [26]. 

This work is partially supported by the convention de cooperation between 
CNRST-Morocco/GRICES-Portugal 681.02/CNR, and by the PROTARS IIP grant 
D16/04. 
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